Singlet-triplet avoided crossings and effective g factor 
versus spatial orientation of spin-orbit-coupled quantum dots 
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We study avoided crossings opened by spin-orbit interaction in the energy spectra of one- and two- 
electron anisotropic quantum dots in perpendicular magnetic field. We And that for simultaneously 
present Rashba and Dresselhaus interactions the width of avoided crossings and the effective g factor 
depend on the dot orientation within (001) crystal plane. The extreme values of these quantities 
are obtained for [110] and [110] orientations of the dot. The width of singlet-triplet avoided crossing 
changes between these two orientations by as much as two orders of magnitude. The discussed 
modulation results from orientation-dependent strength of the Zeeman interaction which tends to 
polarize the spins in the direction of the external magnetic field and thus remove the spin-orbit 
coupling effects. 

PACS numbers: 73.21.La 
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I. INTRODUCTION 

Spin-related phenomena in few-electron quantum dots 
have been under an extensive investigation during the 
past decade. The studies covered spin relaxation in- 
volving spin-orbit interaction and phonon emission, 1 spin 
dephasing 2 due to coupling to nuclear spins 3 as well as 
the spin exchange interaction. 4 Besides the fundamen- 
tal interest these studies were motivated by an idea to 
implement quantum computation^ on spins of separate 
electrons confined in quantum dot arrays. The spin-orbit 
(SO) coupling is considered for spin manipulation within 
the orbital degrees of freedom 6 as well as in the context of 
the anisotropy^ of the exchange interaction for quantum 
gating. 8 

Few-electron systems confined in circular quantum 
dots undergo ground-state angular momentum transi- 
tions in external magnetic field (B)£~— For the electron 
pair these transitions are observed only in presence of 
the electron-electron interaction and are accompanied by 
spin transitions with ground-state changing between sin- 
glet and triplet depending on the parity of angular mo- 
mentum quantum number ^ Singlet-triplet ground- state 
transitions in two-electron quantum dots are observed in 
charging experiments* 1 ^ In elliptical quantum dots as well 
as in double dots^ 3 - the angular momentum transitions 
are replaced by ground-state parity symmetry transfor- 
mations still accompanied by singlet-triplet transitions. 

The SO coupling mixes the eigenstates of opposite par- 
ities and spin. In presence of the SO interaction the 
singlet-triplet transition occurs through an avoided cross- 
ing that for planar quantum dots was discussed in a num- 
ber of recent theoretical papers^ The SO interaction 
usually introduces energetically weak effects hence the 
singlet-triplet avoided crossing is narrow and difficult to 
observe experimentally. The first observations of singlet- 
triplet avoided crossings due to SO coupling were per- 
formed in transport experiments on quantum dots formed 



in gated In As quantum wires ^ A transport experiment 
resolving this avoided crossing in a planar structure was 
performed only recently^ on a single InAs self-assembled 
quantum dot. 

The SO interaction appears due to inversion asymme- 
try of the crystal lattice (Dresselhaus 2 ^ coupling) and/or 
of the nanostructure (Rashba 2 ^ coupling). The resulting 
SO Hamiltonian is not invariant with respect to rota- 
tions within the plane of confinement. The anisotropy of 
SO interaction was investigated by observation of singlet- 
triplet avoided crossing for rotated external magnetic 
field vector J£ The study of Ref. [2 extends the pre- 
vious workii in which the spin splittings were controlled 
by orientation of the external magnetic field superposing 
the effective magnetic fielcU^ introduced by SO coupling. 
It was also demonstrated^ that in presence of the SO 
coupling the non-linear Kondo conductance depends on 
the orientation of external magnetic field. 

In this paper we consider a planar anisotropic quantum 
dot in a perpendicular magnetic field and demonstrate 
that the anisotropy of SO interaction can be used for 
tuning the width of the singlet-triplet avoided crossing 
by spatial orientation of the dot. This tunability appears 
provided that both SO coupling types are present. The 
discussed effect results from dependence of the effective 
strength of the Zeeman interaction on the quantum dot 
orientation within (001) crystal plane. The Zeeman inter- 
action tends to polarize electron spins in the direction of 
the magnetic field. A complete polarization amounts in 
removal of the SO coupling effects. The extent of the spin 
polarization - and thus also the effective Lande factor— 
- vary with the dot orientation. For similar values of 
SO coupling constants the width of the avoided cross- 
ing changes by two orders of magnitude between a few 
fieV to about 0.5 meV. The dependence of the width of 
singlet-triplet avoided crossing on spatial orientation of 
the dot is present for any form of confinement potential 
(quantum well or parabolic profile) for both single and 
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double quantum dots. 



II. THEORY 



A. Hamiltonian 



We consider a quantum dot defined within (001) plane 
with x and y axes identified with [100] and [010] crystal 
directions, respectively. The magnetic field is oriented 
parallel to the growth [001] direction (z). We adopt a 
two-dimensional approximation assuming that the con- 
finement potential is separable into a vertical and planar 
components W(r) = V z (z) + V(x,y) and that the ver- 
tical confinement is much stronger than the planar one. 
Under these conditions the contribution of states excited 
in the vertical directions introduced by the spin-orbit- 
coupling and by the electron-electron interaction is neg- 
ligible. The two-dimensional single-electron Hamiltonian 
takes the form>^ 



H 



hi + ^gfi B Ba z 



I A 



H 



BIAi 



(i) 



where h = (^^r + V(x,y)^ is the spatial Hamiltonian, 

1 is the identity matrix, k = -i(JL, JL) + £(A x ,A y ) , 
g stands for the Lande factor, and Hsia and Hbia de- 
scribe Rashba 22 ( structure inversion asymmetry) and 
Dresselhaus 2 ^ (bulk inversion asymmetry) SO interac- 
tions. We use the symmetric gauge with A x = — , 
A y = x^. The two-dimensional Rashba interaction 
is composed of the diagonal and linear parts Hsia — 



n SIA 



and 
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diag 
SIA 



with 
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a(a x k y - (T y k x ), 



(2) 



777* = 0.0465ra , and the Lande factor g = —8.97. For 
V z in form of a infinite quantum well of height d the 

2 

two-dimensional Dresselhaus constant is f) = 73^ ^ • For 
d = 5.42 nm we have /? = 10.8 meV nm. The two- 
dimensional Rashba constant achieves this value when 
(an external or built-in) vertical electric field acquires 
188.8 kV/cm. 

Below we consider a model confinement potential 



V c (x',y') 



V 



(6) 



where Vb = 50 meV is assumed for the depth of the 
quantum dot. The exponent /i = 10 is applied for which 
the potential profile has a form of a rectangular potential 
well with smoothed boundaries. We take 2K = 40 nm 
as the smaller length of the dot and the larger length is 
taken 2L = 90 nm, unless stated otherwise. The primes 
standing in Eq. (j6j) are referred to the crystal directions 
x and y by 



y 



xcos(<fi) — ys'm((f)) 
xs'm((j)) + 2/cos(0). 



(7) 



The orientation of the dot with respect to the crystal 
directions is displayed in Fig. 1 for </> = 7r/4. The effects 
discussed below remain qualitatively the same for other 
profiles of the dots. At the end of next Section we present 
also results for elliptical parabolic potential as well as for 
a double dot potential. 



B. Method 

The eigenstates of the single-electron Hamiltonian (pQ) 
are calculated in a basis of multicenter Gaussian func- 
tions which is a precise tool for treatment of confinement 
potentials of arbitrary or no symmetry 25 
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The two-dimensional coupling constant a in Eq. ([2]) is re- 
lated to the bulk coupling constant asD as a = a^D (^-)? 
where the average value is calculated for the ground- 
state wave function in the growth direction. The Dres- 
selhaus interaction contains a linear and cubic terms 
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BIA 



= H 
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Hbia = Isd [(Tyk y kl - a x k x k 2 y ] 



(4) 



(5) 



where 73 £> is the bulk coupling constant and the two- 
dimensional constant is defined by /3 = 73£>(& 2 ). We 
adopt the material parameters for an Ino.5Gao.5As quan- 
tum dot, asD = 0.572 nm 2 (see Ref. l23f ) and 73 £> = 
32.2 meVnm 3 (see Ref. [24]), the electron effective mass 



exp 



(r - R fc f 
2a 2 



ieB ( v 



yx k ) 



(8) 

where the centers of Gaussians Rk = (X/^Yfc), are dis- 
tributed on a rectangular array,— and the localization 
parameter a is optimized variationally. In Eq. (|8j) s = ±1 
and Xs are eigenstates of the Pauli matrix <j z . 

The two-electron states are found by the exact diag- 
onalization approach, which uses the basis of the anti- 
symmetrized products of operator (pQ) eigenstates 
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N N 

J2 E (<wi)<M2) 

fl—l V=jJL-\-l 



^(2)^(1)). (9) 



for diagonalization of the two-electron Hamiltonian H2 = 
H(l) + H(2) + 2—^ — (e = 13.55 is taken for the dielec- 

v ' v / 1 47reeori2 v 

trie constant). For 2K = 40 nm and 2L = 90 nm we 
use a basis of 25 x 25 centers, which gives 1250 elements 
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including the spin degree of freedom. In the two-electron 
calculations we take N = 30 lowest-energy single-electron 
spin-orbitals which produces a basis of 435 elements and 
guarantees the convergence of the variational result bet- 
ter than 1/ieV. 
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FIG. 1: The contour shows a quantum dot of width 2K = 
40 nm and length 2L — 200 nm placed along [110] crystal 
direction [</> = ty/4 in Eq. (7)]. With the colors we plotted the 
values of the cosine in the integrand of Eq. (fT3|) for a = 10.8 
meV nm. 
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FIG. 2: (a) The dots show the eigenvalues of Hamiltonian 
(flip and the lines present the results of diagonalization of the 
exact Hamiltonian in function of the external magnetic field 
for a = (3 = 10.8 meV nm, 2K 40 nm and 2L 350 
nm. The results for [110] ([110]) orientation of the dot are 
displayed in red (black) . (b) The spin of the two lowest energy 
levels for both orientations of the dot. The dots in (c) and 
(d) show the results of diagonalization of four by four matrix 
Hamiltonian with basis including the first excited state (see 
text). 



III. RESULTS AND DISCUSSION 

A. Effective g factor and orientation of the dot 

In order to explain the dependence of the strength of 
the Zeeman interaction on the orientation of the dot - 
which underlies the results to be presented below - let 
us consider the special case of equal linear SO coupling 
constants a = j3 = 10.8 meV nm and neglect the cu- 
bic Dresselhaus and diagonal Rashba terms of SO in- 
teraction, which are small anyway^ We consider the 
approximate Hamiltonian for B = defined as Ho = 
hi + H l gf A + H l Bj A . Hq commutes with the operator 
of [110] spin component and SO coupling shifts the en- 
tire electron energy spectrum by a constant quantity 27 
En = £n — , where En and cn denote energy 

eigenvalues with and without SO coupling, respectively. 
For B = the SO coupled eigenfunctions of H opera- 
tor ^/n± corresponding to ±h/2 spin eigenvalues in the 
[110] direction are related to orbital eigenfunctions ipx 
that are obtained in the absence of SO coupling as 

<l>N±(x,y) = ( ±e -W4 ) VN(x,y)e^-^( x+y) . 

(10) 

The magnetic field vector parallel to the growth direc- 
tion introduces the Zeeman interaction with cr z matrix 
to the Hamiltonian, and the [110] spin component is 



no longer a good quantum number. Lets try to di- 
agonalize the Hamiltonian including the Zeeman effect 
H z = Hq + ^giiBB(7 z taking Hq eigenstates ([TO]) as the 
basis. The shortest reasonable basis contains two de- 
generate ground- state wave functions (j)\± corresponding 
to opposite spin orientations and the same orbital wave 
function ip\ . The matrix of H z operator takes the form 



H z 



(ii) 



where both diagonal terms are (0i±|#b|0i±) = E\ 

9 2 

a h2 m , and the off-diagonal ones are 



(cj> 1± \H B \cj> lT ) = 

y f i B Bj\ ( p 1 (x,y)\ 2 e( ±i "^ 1{x+y) )dxdy. (12) 

For potentials with an in-plane inversional symmetry 
that are considered in this paper the matrix element (fT2j) 
is real and given by 



J \cpi(x, y)\ 2 cos 



W2 



(x + y) dxdy. (13) 



Figure 1 shows the plot of the cosine term in the inte- 
grand. The argument of the cosine has a fixed orienta- 
tion with respect to the crystal directions and changes 



sign along [110] direction with a period of Xso = 2am* 
(for the applied parameters Xso = 238.4 nm). For the 
quantum dot of length L = 200 nm oriented along [110] 
(see Fig. 1) the cosine has the same sign within the 
quantum dot area. For this orientation the off-diagonal 
terms of Hamiltonian matrix (11) are the largest. On 
the other hand for the dot oriented along [110] the sign 
of the integrand oscillates within the quantum dot area 
and the off-diagonal terms are necessarily smaller. The 
off-diagonal matrix elements mix the cr x - y eigenstates 
leading to alignment of the spin along the direction of 
the field (z). Therefore, the spin polarization due to the 
Zeeman effect should be the strongest for [110] dot ori- 
entation and the weakest for the dot oriented along [110] 
crystal direction. 

The results of diagonalization of matrix Hamiltonian 
(fTTj) is displayed in Fig. 2(a) for the dot oriented along 
[110] (black dots) and [110] (red dots) directions. In Fig. 
2(b) we displayed the average spin for the two lowest- 
energy states for both orientations of the dot. The lines 
in Fig. 2 show the results of the exact diagonalization 
with basis given by ((8]). The eigenvalues of matrix ([TT]) 
quite well reproduce the exact energy levels and the av- 
erage spin.— According to the intuition given by Fig. 1 
the electron spin reacts to the application of the external 
magnetic field in a more pronounced manner for [110] ori- 
entation of the dot than for the perpendicular orientation 
[110] . As the non-zero magnetic field lifts the ground- 
state degeneracy the ground state (the first excited state) 
becomes nearly spin-up (spin-down) polarized. Polariza- 
tion of the spin by infinitesimal B for the dot oriented 
along [110] direction is much weaker and increases for 
higher fields. This increase [black lines in Fig. 2(b)] is 
not very well reproduced by the two-element basis (black 
dots). Inclusion of the first excited scalar wave function 
cf2 to the approximate calculation gives four basis ele- 
ments of type (10). The results are displayed in Figs. 
2(c) and 2(d). The four-element basis reproduces also 
the excited energy levels and an improvement of the de- 
scription of s z (B) dependence is obtained particularly for 
the [110] orientation. 

The extent of the spin polarization that varies with the 
dot orientation results in the dependence of the Kramers 
multiplet splitting induced by weak magnetic fields. This 
in turn leads to orientation dependence of the effective 
g* factors^ that in the experiments are estimated by the 
splitting of energy levels by weak magnetic field. We 
estimate the effective factor by 



lim 



E2 — E\ 



B-,0 fi B B 



(14) 



which for the antidiagonal [110] orientation of the dot 
gives g* = —8.7 (quite close to g = —8.97) and for the 
diagonal orientation [110] only g* = —2.5. 



B. Single-electron results 

Let us now consider the results obtained by diagonal- 
ization of Hamiltonian (pQ) with basis (|8j). Figure [3] shows 
the energy levels and mean spin z components for var- 
ious length of the dot 2L for the diagonal [110] (black 
color) and the antidiagonal [110] (red color) dot orien- 
tation. For the dot which is close to the square profile 
[2K = 40 nm and 2L = 50 nm, see Fig. 02(a)] we plot- 
ted the results without SO coupling by the green dotted 
line (shifted down on the energy scale by 0.142 meV). 
The first and second energy levels correspond to oppo- 
site parity and spin. For B below the avoided crossing 
the first excited state is of the even parity with spin ori- 
ented down and the second excited state is of the odd 
parity with spin oriented up. The SO coupling opens an 
avoided crossing between these two energy levels. This 
avoided crossing is wider for the diagonal and thinner for 
the antidiagonal dot orientation [Fig. E3[a)]. The width of 
this avoided crossing is determined by an extent to which 
the SO coupling entangles the spin and orbital wave func- 
tions. The width varies more strongly with the dot ori- 
entation when the anisotropy of the dot is enhanced, i.e. 
for larger length of the dot [Fig. 02(a,c,e,g)], particularly 
when it becomes comparable to Xso- For the antidi- 
agonal orientation of the dots the strong Zeeman effect 
quickly polarizes the electron spin and thus removes the 
SO coupling effects from the energy spectrum. The en- 
ergy spectra obtained without SO coupling are close to 
the ones obtained for the antidiagonal orientation of the 
dot - see Fig. 02(a) and (e). 

Figure [4] shows the dependence of the effective g* factor 
on the orientation of the dots for 2L = 90 nm and 200 
nm. The g* factor acquires maximal (minimal) absolute 
value for the antidiagonal (diagonal) dot orientation. For 
non-equal coupling constants variation of g* is reduced, 
and disappears for a single type of SO coupling present. 
Note that for the antidiagonal orientation of the dot the 
same g* is obtained for both L considered. 

For completeness in Fig. 03(b) we displayed the mean 
value of the z component of the spin in function of the dot 
orientation angle for a = /?, 2L = 90 nm and two values 
of the magnetic field B = 1 T (red) and B = 3 T (black 
curves). The spin polarization is the largest for the <p = 
7r/4 and the smallest for <p = — 7r/4. Fig. 03(a) shows the 
probability density of spin-up and spin-down components 
for the first excited state. Both components of the wave 
function have opposite spatial parities. For <fi = 7r/4 the 
spin-up component is not visible in the applied contour 
scale. 

The spatial orientation of the dot has a significant in- 
fluence on the SO related avoided crossings only when 
both types of the coupling are of comparable strength. 
Fig. [6] shows the results for dominant Dresselhaus term 
a = 0.2/3 for 2L = 90 nm. The dependence of the 
width of the avoided crossing on orientation is qualita- 
tively the same as for a = f3 [see Fig. 03c, d)], only much 
weaker. For comparison the energy spectrum for a = 
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2L=50 nm 2L=90 nm 2L=200 nm 2L=350 nm 




B[T] B[T] B[T] B[T] 

FIG. 3: Energy spectrum (upper row) and mean values of the z component of the spin (lower row of plots) for a single electron 
in a quantum dot of width 2K and various lengths 2L. The results for the dot oriented along [110] and [110] are given in black 
and red, respectively. In (c) and (d) we additionally display the results for the dot placed along the [010] direction. Equal 
Rashba and Dresselhaus linear coupling constants were assumed a — f3 — 10.8 meV nm. In (a) and (e) the green dotted lines 
show the energy spectrum in the absence of SO interaction shifted down on the energy scale by 0.142 meV. 
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FIG. 4: Effective g* factor estimated by Eq. (14) in function 
of the spatial orientation of the dot for 2K = 40 nm, 2L = 90 
nm (solid lines) and 2L = 200 nm (dashed lines) for f3 = 10.8 
meV nm and various values of a. 



and j3 = 10.8 meV nm is given in Fig. [6ja) with the black 
dotted line. For a single type of SO coupling present the 
same energy spectrum is obtained for any dot orientation. 



C. Two-electron results 

The two-electron spectrum without SO coupling is 
given in Fig. [7Ja) by the green dotted lines. When the 
magnetic field is swept across B = 1.1 T the ground state 
changes from the singlet with even spatial parity to the 



FIG. 5: Results for 2L = 90 nm and a = f3 = 10.8 meV 
nm. (a) Spin densities for the first excited-state calculated 
for subsequent (j) values (— 7r/4, 0, ty/4 and ty/2) for B = 3T. 
(b) (s z ) for B = IT (B = 3T) plotted in function of angle 
with black (red) lines respectively. 



spin-up polarized triplet with odd spatial parity. The 
dashed and solid lines in Fig. [71 correspond to SO cou- 
pling for a = f3 = 10.8 meV nm. The results for [110] 
and [110] dot orientations are plotted with the black and 
the red lines, respectively. For [110] dot orientation the 
singlet-triplet transition produces a very narrow avoided 
crossing of width 6fieV as compared to the pronounced 
(0.37 meV wide) avoided crossing obtained for [110] ori- 
entation. For [110] orientation the mean values of the 
spin vary smoothly as functions of B [see Fig. [7(b)], 
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2 3 
B[T] 



FIG. 6: (a) Energy spectrum for a single-electron dot with a 
larger length aligned with [010] (blue), [110] (black), and [lTO] 
(red) crystal directions given by solid and dashed lines. The 
dot size is 2K = 40 nm and 2L = 90 nm, f3 = 10.8 meV nm 
and a = 0.2/3. The dotted curve shows the results without the 
linear Rashba term (a = 0). (b) Mean value of s z operator 
for three lowest energy states respectively for [110] and [110] 
orientations. The results for this geometrical parameters and 
a = P were given in Fig. G2c,d). 



while for [110] the mean spin is nearly a bi- valued func- 
tion of magnetic field [see Fig. [7(b)] , which indicates a 
removal of the SO coupling effects. The corresponding 
energy spectrum [red lines in Fig. [7(a)] is very close to 
the spectrum obtained without SO coupling (green dot- 
ted lines) up to a constant energy shift [the energy levels 
without SO coupling are shifted down by 0.285 meV in 
Fig. [7(a)]. Similar fact was presented above for the single 
electron in Fig. 3(a). 

In Fig. [8(a) we present two-lowest energy levels as 
functions of the orientation of the dot for three values 
of magnetic field: 1/iT (residual £?), 0.5 T - before the 
singlet-triplet avoided crossing - and for 1.1 T - at the 
center of avoided crossing. The average spins are dis- 
played in Fig. [8(b). For the residual magnetic field the 
energy spectrum is independent of the dot orientation 
and the ground-state spin is zero. Nevertheless, a depen- 
dence of the spin of the excited state (three-fold degen- 
erate at B = 0) on the orientation of the dot is notice- 
able. For 0.5 T the energy levels weakly depend on the 
dot orientation, but dependence of the spins is strong. 
The situation is opposite for 1.1 T. In both the cases for 
(j) = 7r/4 the spins approach closest and % - values that 
are found in the absence of SO interaction. 



u L i j 

FIG. 7: (a) Two-electron energy spectrum for single- 
elongated-dot aligned along [110] (black), [110] (red). Green 
dotted curve shows the results without SO coupling shifted 
down by 0.285 meV. (b) Mean value of s z for two lowest en- 
ergy states (solid, short dashed respectively) for [110] (black) 
and [110] (red) dot orientation. The value of a = f3 = 10.8 
meV nm is assumed, 2K = 40 nm and 2L = 90 nm. 
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FIG. 8: (a) Two lowest two-electron energy levels plotted in 
blue for B = 1/iT, in black for B = 0.5T, and in red for 
B = LIT. (b) The s z mean value. Same parameters as in 
Fig. [71 were adopted. 



D. Parabolic and double dots 

The profiles of the confinement potentials depend 
strongly on the type of quantum dots, their size and 
growth conditions. The effects discussed above occur in 
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the low-energy part of the spectrum and appear as func- 
tions of the dot orientation. In order to demonstrate that 
they are not specific to any profile of confinement poten- 
tial we considered also a parabolic quantum dot with an 
elliptical shape and a double dot. 
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FIG. 9: One (a) and two-electron (b) energy spectrum for 
elliptical dot aligned along [110] (black) and [110] (red) for 
the elliptic parabolic confinement potential given by Eq. {15)) . 
Results are obtained for a — f3 = 10.8 meV nm. The insets 
present the equipotential lines. 



N = 2 




FIG. 
(16). 



10: Same as Fig. [9] only for the double dot potential 



Confinement potential of electrostatic quantum dots 29 
is generally parabolic close to its minimum, although 
potential profiles closer to a quantum well can also be 
realized^ We considered the potential in form 



V P (x',y') = 



TOW, 



x x' 2 



tow; 



V 2 , 



(15) 



with %uu x = 2 meV and Tiuu y = 5 meV. The results for x' 
identified with [110] and [110] directions are given in Fig. 
[9J The avoided crossing between the first and second 
excited states of the single-electron [Fig. [9ja)] as well 
as the singlet-triplet [Fig. l2fb)] avoided crossing vary 
strongly with the orientation of the dots in consistence 
with the above discussion for the confinement potential 
given by Eq. ((6]). 

We model a double dot by introducing a barrier in the 
center of the quantum dot 



V d (x',y') = V c (x',y') 



(16) 



(i+[&r) (*+[£]) 



where V c is defined by Eq. (j6j) with 2K = 40 nm and 
2L = 90 nm, V& = 10 meV, and 2B = 10 nm is taken 
for the barrier width. For comparison the results for 
the single dot were presented in Fig. 3(c) for the sin- 
gle electron and in Fig. 6(a) for the electron pair. Both 
the single- and two-electron avoided crossings that were 
discussed above involved mixing of even and odd spa- 
tial parity states by SO interaction. For the double dot 
these states correspond to bonding and antibonding or- 
bit als, respectively. The avoided crossings observed for 
the double dot are considerable thinner than in the single- 
dot case which is due to the introduction of the interdot 
barrier (in the limit of an impenetrable interdot barrier 
the bonding and antibonding orbitals are degenerate). In 
Fig. [10] we find the dependence of the width of avoided 
crossing on the orientation of the dot that agrees with 
the precedent results. 



E. Discussion 

The presented results indicate that the width of SO- 
related avoided crossings can be designed by specific ori- 
entation of the dot with respect to the crystal axes. The 
choice of the orientation has to be done at the sample fab- 
rication stage. In gated quantum dots with confinement 
potential of the electrostatic origin 2 ^ the orientation of 
the dot can be chosen quite arbitrarily by the shape of 
electrodes defined on the sample surface. Orientation 
of quantum dots with structural confinement can also 
be intentionally controlled. For instance, InGaAs/GaAs 
double quantum dots are formed on pre-patterned sub- 
strates along either [110] or [110] directions, 31 for which 
the width of the SO-related avoided crossing acquires ex- 
tremal values. In electrostatic quantum dots with a mul- 
tielectrode setup,— rotation of the confinement potential 
should be possible to realize by voltages applied to the 
gates on a single sample. 

We find that the orientation of the dot influences the 
width of avoided crossings provided that both the lin- 
ear SO coupling constants are similar. The Rashba in- 
teraction constant can be adjusted by external electric 
fields, 32 to match the Dresselhaus constant in particular 
for observation 3 - 3 - of persistent spin helix states in quan- 
tum wells. 34 When only a single type of SO coupling is 
present the orientation of the dots has no influence on 
the energy spectrum. In presence of the SO coupling 
the orbital angular momentum is not a good quantum 
number even for circular confinement potentials. Never- 
theless, both the spin components of single-electron wave 
function do possess a definite albeit different angular mo- 
menta. In consequence the charge density reproduces the 
circular symmetry of the confinement potential^ In po- 
tentials of circular symmetry the charge density becomes 
anisotropic only when both SO coupling types are present 
and additionally the Zeeman effect is introduced by ex- 
ternal magnetic field.— 

The effect of the dot orientation on the energy spec- 
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trum is obtained in external magnetic field. For 5 = 0, 
the dot orientation has no influence on the Kramers- 
degenerate energy spectrum even when both SO coupling 
types are present. For a = /3 the effective magnetic 
field introduced by SO coupling is oriented along the 
[110] direction independent of the orientation of the dot, 
nevertheless its strength is dot-orientation-dependent ^ 
The electron spins precess in the effective magnetic field. 
Therefore, the orientation of the dot does matter for 
the spin-manipulation at zero magnetic field j 27 i 36 even 
though no effect on the energy spectrum is present. 

IV. SUMMARY AND CONCLUSIONS 

We studied the avoided crossings opened by SO in- 
teraction in the single- and two-electron planar (001) 
quantum dots as functions of the external perpendicu- 
lar magnetic field. We demonstrated that the width of 
these avoided crossings can be tuned within a range of 
two orders of magnitude by orientation of the quantum 
dot with respect to the crystal directions. 



The tunability is achieved provided that i) both 
Rashba and Dresselhaus interactions are present with 
comparable values of linear coupling constants (a ~ 0) 
ii) the dot is anisotropic and iii) its larger length is com- 
parable to Xso — 7rh 2 /2am*. The dependence of the 
width of avoided crossings on the orientation of the dot 
results from a different strength of the Zeeman interac- 
tion which more or less efficiently polarizes the electron 
spin. The spin polarization removes the spin-orbital en- 
tanglement from wave functions along with the SO cou- 
pling effects from the energy spectrum. Thus, the dot 
orientation affects simultaneously the width of avoided 
crossings and the effective Lande factor g*. As a gen- 
eral rule, the dot orientations producing large \g*\ values 
correspond to narrow avoided crossings. 
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